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Abstract. A method for the calculation of the luminosity for the proton-nucleus collisions based 
on the quasi-free proton-proton scattering is presented. As an example of application the integrated 
luminosity for the scattering of protons off the deuteron target is determined for the experiment of 
the quasi-free pn pnr] reaction performed by means of the COSY- 1 1 facility. 
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INTRODUCTION 

Due to the lack of the pure neutron targets a meson production in the proton-neutron 
reactions is usually realised e.g. via the proton scattering off the deuteron. This method 
takes advantage of the fact that the neutron's binding energy inside the deuteron is 
negligible relative to the kinetic energy of proton beams needed to produce the meson in 
the nucleon-nucleon collisions In the analysis of such reactions it is assumed that 

the proton is acting only as a spectator and that it is on its mass shell at the moment of 
the collision. 

In this contribution we present the evaluation of the luminosity in such kind of 
experiments by means of the measurement of the quasi-free proton-proton scattering. 
This method permits to take into account automatically the shadowing effects of the 
spectator nucleon, and allows for the determination of the luminosity with a relatively 
small normalization uncertainty thanks to the availability of the precise cross sections 
for the proton-proton elastic scattering determined by the EDDA collaboration [3] up to 
the beam momentum of 3.5 GeV/c. 



LUMINOSITY 

In order to determine the luminosity for the quasi-free proton-neutron reaction we mea- 
sured quasi-free pp ^ pp reaction by detecting in coincidence both scattered protons. 
Here we will describe the measurement performed by means of the COSY- 11 appara- 
tus [0, S 0, 0] shown schematically in Figure [H (left). The recoil proton gives signal 
in the scintillator detector S4 and subsequently reaches the granulated silicon detector 
^^mon^ while the forward scattered proton is registered by the stack of drift chambers 
Dl and D2 and scintillator array SI. For triggering of the elastic scattered events the 
coincidence between signals from the SI and S4 scintillators was required. 




FIGURE 1. (left) Schematic view of the COSY- 11 detection setup. Only the detectors used for the 
measurements of the quasi-free proton-proton scattering are shown. Dl and D2 denotes the stack of two 
drift chambers. Simon is the granulated silicon detector. SI and S4 stand for the scintillator detectors, (right) 
Relation between the center-of-mass scattering angle and the position in the S 1 detector for the quasi-free 
proton-proton scattering at Pbeam = 2.075 GeV/c as has been obtained in the Monte-Carlo simulations. 



In the case of free proton-proton scattering the luminosity could be determined as a 
normalization constant between the measured angular distribution of the cross section 
and the corresponding spectrum known from the previous experiments. For the quasi- 
free reaction the evaluation becomes more complicated due to the Fermi motion of 
nucleons inside the nucleus (Fig. [2] (left)). Since the direction and momentum of the 
bound nucleon may vary from event to event this implies that the direction of the 
center-of-mass velocity of the colliding nucleons as well as the total available energy 
for the reaction may also vary from event to event. Therefore, protons registered in the 
laboratory under a given scattering angle or at a given part of the detection system, 
correspond to the finite range of scattering angles in the proton-proton center-of-mass 
frame (see Fig.[I](right)). This implies that the experimental angular distributions cannot 
be directly compared to the literature values, and instead an evaluation of the luminosity 
requires simulations taking into account the Fermi motion of the nucleons, and the 
variations of differential cross sections for the elastic scattering as a function of the 
scattering angle and energy. 

For each simulated event we know the generated Fermi momentum of the nucleon, as 
well as the scattering angle of protons in their center-of-mass system. This permits us 
to assign to each event a weight corresponding to the differential cross section, which is 
uniquely determined by the scattering angle and the total collision energy s. 

For the intuitive illustration of the size of the momentum spread caused by the 
Fermi motion instead of the total center-of-mass energy we may equivalently consider 
the effective beam momentum as it is seen from the nucleon inside the nucleus. The 
distribution of the effective beam momentum depends on the value of the proton beam 
momentum and as an example in Figure [2l (middle) we present it for the value of 
2075 MeV/c as which used for the measurement of the quasi free pn pni] reaction [@, 
@]. It is important to note that the distribution of the equivalent beam momentum ranges 
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FIGURE 2. (left panel) Nucleon momentum distribution inside the deuteron according to the Paris [[K 
(full line) and CDBONN |11] (dotted line) potentials, (middle panel) Distribution of the beam mo- 
mentum as seen by the proton bounded inside a deuteron hit by the beam proton with momentum of 
Pbeam = 2.075 GcV/c. (right) Bilinear interpolation of the differential cross section ^{pteam^ ^*)- 



from about 1.5 GeV/c up to circa 2.5 GeV/c. In this momentum range the cross sections 
for the proton-proton elastic scattering vary significantly [3.] and therefore this effect 
cannot be neglected. 

In the following we will be more specific and will describe the derivation of the 
luminosity in terms of formulae exploited in the analysis. For a free proton-proton 
scattering we could measure the number of events - AA/^(0,0) scattered into the solid 
angle AQ.{6,^) around the polar and azimuthal angles 6 and respectively. In this case 
the angles in laboratory and in the center-of-mass systems are univocally related to each 
other. With the known differential cross section [3] for proton-proton scattering into that 
particular solid angle, and having known the value of the solid angle AQ.{6,^) from the 
Monte-Carlo simulations the luminosity can be calculated according to the formula: 

L= ^'^/> . ,1) 

An(e,,/i)^(9,*) 

In the case of quasi-free proton-proton scattering the number of elastically scat- 
tered protons AN into a solid angle AQ.{6iab, ^lab) is proportional to L - the integrated 
luminosity over the time of measurement, and also to the inner product of the dif- 
ferential cross section for scattering into the solid angle around 0* and </)* angles - 
^(0*,(/)*,;7/7,0/7,(/)/7) - and the probability density of the distribution of the Fermi mo- 
mentum /(/7/7, Of^(^p): 

ANejcp{Aa{eiab,(^lab)) = 

JAaiOiabAab) 

The angles 0* and (/)* are expressed in the proton-proton center-of-mass system, 
while the angles Oiab and (^lab are considered in the laboratory system. In the case of 
the complex detection geometry with a magnetic field a solid angle corresponding to a 
given part of the detector cannot in general be expressed in a closed analytical form. 



Therefore, integral in Equation [2] must be computed using the Monte-Carlo simulation 
programme, containing the exact geometry of the detection system and taking into 
account bending of particles trajectories in the magnetic field as well as detection and 
reconstruction efficiencies. For the evaluation of a given event by the Monte-Carlo 
programme first we choose randomly a momentum of a nucleon inside a deuteron 
according to the Fermi distribution [10] (Fig. [21 (left)). Next the total energy s for the 
proton-proton scattering and the vector of the center-of-mass velocity are determined. 
Then, we generate isotropically a momentum of protons in the proton-proton center- 
of-mass frame. Further on, according to the generated angle and the total-energy s (or 
equivalently an effective beam momentum seen by the struck nucleon) we assign to the 
event a probability equal to the differential cross section. Next, the momenta of protons 
are transformed to the laboratory frame and are used as an input in the simulation of the 
detectors signals with the use of the GEANT computing package. 

The differential cross sections ^(Q*? with s being dependent on pp^Op,(^p, 
and Pbeam^ ^nd with 0* and (/)* denoting the scattering angles in the proton-proton 
center-of-mass system were calculated using the cross section data base for pp ^ 
pp reaction ^ H For this purpose we have applied a bilinear interpolation in the 
momentum-scattering angle plane and calculated differential cross section according 
to the formula: 

tU^{pleam^ S*'^) + (1 - t)u^{pl^^, 0*'^), (3) 

where variables t and u are defined in Figure [21 (right). These cross sections were used 
as the weights of the elastically scattered events in the Monte-Carlo calculations. 

/^exp from Equation [21 can be determined as a number of elastically scattered protons 
registered in a given part of the detector system. In order to calculate the integral on 
the right hand side of this equation we simulated A/q events according to the procedure 
described above. Due to the weights assigned to the events the integral is not dimension- 
less and its units correspond to the units of the cross sections used for the calculations. 
The number obtained from the Monte-Carlo simulations must be then normalized such 
that the integral over the full solid angle equals to the total cross section for the elastic 
scattering averaged over the distribution of the total reaction energy s resulting from the 
Fermi distribution of the target nucleon. In the absence of the Fermi motion it should 
be simply equal to a total elastic cross section for a given beam momentum. This means 
that we need to divide the resultant integral by the number of generated events A/q and 
multiply it by the factor of Ik, A factor 27i comes from the normalization of the differ- 
ential cross section ^{pbeam^ 0*, (/)*), regarding the fact that protons taking part in the 
scattering are indistinguishable. Hence, the formula for the calculation of the integrated 



^ EDDA group has gathered the data for the excitation functions ^{0^^ ^Pbeam) for the elastic pp ^ 
pp process at 108 different proton kinetic energies, ranging from 240 MeV up to 2577 MeV. In the 
measurements the center-of-mass scattering angles of protons (0*) from 30° up to 90° have been covered. 
Both the kinetic and angular ranges are sufficient to cover our needs in calculating the corresponding 
differential cross sections. 



luminosity for the quasi-free reaction reads: 

(4) 

where the normahzation constant No/2n is subject to the Monte-Carlo method used for 
the integral computation. 

EXAMPLE OF APPLICATION 

We have applied the above described method for the evaluation of the luminosity for the 
measurement of the pn pni] reaction with the COSY- 11 facility [f8l.[9l.[l3[]. 

Events corresponding to the elastically scattered protons have been identified on 
the basis of the momentum distributions. The momentum of the fast scattered proton, 
whose trajectory has been registered in the drift chambers, can be reconstructed, and 
the transversal versus the parallel momentum component may be plotted as it is done in 
left panel of Figure [3l The signal from the elastic scattered protons appears as an clear 
enhancement around the expected kinematical ellipse. 




FIGURE 3. Parallel versus transversal momentum component of the reconstructed fast proton momen- 
tum as obtained in the experiment (left) and in the simulations (right). Superimposed lines correspond to 
the expected kinematical ellipses. 

For the calculations of the integrated luminosity, the part of the SI detector available 
for the elastically scattered protons was divided into four subranges. In order to separate 
the background from the multi particle reactions for each subrange, the distribution of 
the distance of the points to the kinematical ellipse was determined. The result obtained 
for different sections of the SI detector are presented in Figure IH 

A linear background cut El has been performed and subsequently the true scattering 
yields into a given range of the SI detector have been calculated. 

For the determination of the integral of Equation [21 an Nq = 10^ quasi free pp ^ pp 
events have been simulated, and the response of the detectors has been generated using 



^ The main source of the background are the accidental coincidences originating from the production and 
scattering processess. 
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FIGURE 4. Projection along the expected kinematical ellipse of the experimental event distribution 
from Figure [3] (left) for four subranges of the SI detector. 




distance to ellipse [ a.u. ] distance to ellipse [ a.u. ] distance to ellipse [ a.u. ] distance to ellipse [ a.u. ] 

FIGURE 5. Projection along the expected kinematical ellipse of the simulated event distribution from 
Figure [3] (right) for four subranges of the SI detector. 



the GEANT-3 code based simulation programme, maintaining the experimental condi- 
tions of beam and target [14]. Subsequently, the simulated events have been analyzed in 
the same way as the experimental data. Figure [5] shows the simulated spectra analogous 
to the experimental distributions of Figure |4l As described in the previous sections each 
entry in the shown histograms was weighted according to the differential cross sections 
and hence the integral of these spectra normalized to number of simulated events and 
multiplied by a factor of 2n can be substituted for an integral in Equation [2l Thus the in- 
tegrals of experimental (Fig.© and simulated histograms (Fig. [5]) applied in Equation [21 
permitted to determine the luminosity for each of the subrange of SI detector separately. 
The weighted average over the four quoted values of the integrated luminosity equals to 
L = (2.08±0.03)-10^^ cm-2. 

One source of the systematic error may be attached to the background misidentifica- 
tion. Since the line describing the background can be well defined on both: the left and 
the right side of the scattering peaks, we assumed, conservatively, that the systematic 
error due to the assumption of the linearity of the background in the way as presented in 
Figure |4] is less than 20%. The background events constitute around 15% of all events, 
hence the overall systematic error due to the background subtraction is not greater than 
3%. 

Another source of the systematic error originates from the assumption of the bilinear 
approximation of the cross section shown in Figure [21 (right). To estimate this systematic 
uncertainty we made the zeroth-order assumption in which instead of the interpolation 
we took the cross section value from the closest data point in the effective proton beam 
momentum-scattering angle plane. Higher order approximations of the differential cross 



section given by Equation [3] should be not greater than the difference between the zeroth 
order approximation explained above and the bilinear approximation. The performed 
calculations shows that this difference is smaller than 0.2%. 

Taking into account the normalisation error of the EDDA differential cross sections 
(equal to circa 4% [|3|, ll5]), the systematical error originating from the assumption of 
the potential model of the nucleon bound inside the deuteron (equal to about 2%), and 
the two abovementioned sources of the systematical errors, we estimated the overall 
systematical error of the integrated luminosity to be not greater than 9.2%. 
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